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Abstract 

We review the results from Mariner 10 regarding Mercury's gravity field and the results from radar 
OO . ranging regarding topography. We discuss the implications of improving these results, including a 

I determination of the polar component, as well as the opportunity to perform relativistic gravity tests 

with a future Mercury Orbiter. With a spacecraft placed in orbit with periherm at 400 km altitude, 
. apherm at 16,800 km, period 13.45 hr and latitude of periherm at +30 deg, one can expect a significant 

I improvement in our knowledge of Mercury's gravity field and geophysical properties. The 2000 Plus 

■ mission that evolved during the European Space Agency (ESA) Mercury Orbiter assessment study 
(Hechler, 1994) can provide a global gravity field complete through the 25th degree and order in spherical 

O I harmonics. If after completion of the main mission, the periherm could be lowered to 200 km altitude, 

' the gravity field could be extended to 50th degree and order. We discuss the possibility that a search 

■ for a Hermean ionosphere could be performed during the mission phases featuring Earth occultations. 
Because of its relatively large eccentricity and close proximity to the Sun, Mercury's orbital motion 

k>( ^ provides one of the best solar-system tests of general relativity. Consequently, we emphasize the number 

j_j ■ of feasible relativistic gravity tests that can be performed within the context of the parameterized post- 

^ Newtonian formalism - a useful framework for testing modern gravitational theories. We pointed out 

that current results on relativistic precession of Mercury's perihelion are uncertain by 0.5 %, and we 
discuss the expected improvement using Mercury Orbiter. We discuss the importance of Mercury Orbiter 
for setting limits on a possible time variation in the gravitational constant G as measured in atomic 
units. Moreover, we mention that by including a space-borne ultrastable crystal oscillator (USD) or 
an atomic clock in the Mercury Orbiter payload, a new test of the solar gravitational redshift would 
be possible to an accuracy of one part in lO'^ with a USO, and to an accuracy of one part in 10^ with 
an atomic standard. With an atomic clock and additional hardware for a multi-link Doppler system, 
including Doppler extraction on the spacecraft, the effect of Mercury's gravity field on USD's frequency 
could be measured with an accuracy of one part in 10^. We discuss other relativistic effects including 
the geodetic precession of the orbiter's orbital plane about Mercury, a planetary test of the Equivalence 
Principle (Nordtvedt effect), and a solar conjunction experiment to measure the relativistic time delay 
(Shapiro effect). 

Key Words: Mercury, Mercury Orbiter, Radio Science, Mercury gravity. Mercury ionosphere, Rela- 
tivistic Gravitation 
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1 Introduction 



Mercury is the least explored of the terrestrial planets. Its known global geophysical properties include 
its mass, radius, and gravitational coefficients J2 and C22 (Anderson et al, 1987). Recently, ground- 
based radar ranging data and two range determinations from Mariner 10 have been used to determine 
both the magnitude and direction of the equatorial offset of the center of Mercury's figure from its 
center of mass, and to determine the orientation and major and minor axes of the ellipse which best 
fits Mercury's equatorial shape (Anderson et al, 1995). Displacements of the center of figure from 
the center of mass are known for Earth (Balmino et al., 1973), Moon (Kaula et al., 1972; Bills and 
Ferrari, 1977; Smith et al, 1995), Venus (Bindschadler et al, 1994), and Mars (Standish, 1973; Bills 
and Ferrari, 1978). The global equatorial topography of Mercury from ground-based radar ranging is 
available (Harmon et al, 1986; Harmon and Campbell, 1988; Pitjeva, 1993). The known shape and 
orientation of Mercury's equatorial figure and the displacement of this figure from the planet's center of 
mass place important constraints on the structure of Mercury's interior. The magnitude of the center 
of figure offset implies an excess crustal thickness of 12 km or less, comparable to the Moon's excess. 
By comparing the equatorial ellipticity with the Mariner 10 gravity coefficient C22, and assuming Airy 
isostatic compensation, Anderson et al. (1995) conclude that Mercury's crustal thickness is in the range 
100 to 250 km. A similar calculation for the Moon, using C22 = 2.2 x 10~^ and (o — b)/a = 1.4 x 10~^, 
gives a mean crustal thickness of about ~ 72 km (Schubert et al., 1995). Because the radar ranging data 
are restricted between 12° and - 10° in latitude, Mercury Orbiter can provide additional information on 
geophysical parameters aligned with the polar axis. It can also provide a global gravity field complete 
through the 25th degree and order in spherical harmonics (A polar orbit with periherm at 400 km 
altitude, apherm at 16,800 km, period = 13.45 hr, latitude of periherm at + 30 deg). If after completion 
of the main mission, the periherm could be lowered to 200 km altitude, the gravity field could be 
extended to 50th degree and order. 

Radio propagation measurements using spacecraft signals have also provided scientific information on 
the atmospheres, ionospheres, rings, tori, and surfaces of essentially every planetary system encountered 
in the brief history of solar system exploration (e.g. Tyler, 1987). A search for a Hermean ionosphere 
could be performed with a Mercury Orbiter during the mission phases featuring Earth occultations. 

Because of its relatively large eccentricity and close proximity to the Sun, Mercury's orbital mo- 
tion provides one of the best solar-system tests of general relativity. Its excess perihelion precession, 
amounting to 42.98 arcscconds per century, is one of the three classical tests of the theory. Coherent 
Doppler and ranging measurements using a Mercury Orbiter could provide improved orbital information 
and hence, more sensitive tests. A number of predictions of alternative gravitational theories might be 
detected as well, or failing a detection, tighter constraints could be placed on permissible deviations 
from general relativity. 

By including a space-borne ultrastable crystal oscillator (USO) or an atomic clock in the Mercury 
Orbiter payload, a new test of the gravitational redshift would be possible. The effect of solar gravity 
on the oscillator's frequency could be measured to an accuracy of one part in 10^ with a USO, and to 
an accuracy of one part in 10^ with an atomic standard. With an atomic clock and additional hardware 
for a multi-link Doppler system, including Doppler extraction on the spacecraft, the effect of Mercury's 
gravity field could be measured with an accuracy of one part in 10^. 

This paper summarizes and discusses the different radio science experiments proposed for a future 
Mercury Orbiter mission. The outline of the paper is as follows: Section II discusses a present knowledge 
of Mercury's gravity field and future improvements offered by Mercury Orbiter. In section III, we discuss 
the possibility of a search for a Hermean ionosphere. Mercury ranging data are reviewed in section IV. In 
section V, we discuss feasible relativistic gravitational tests using the satellite in orbit around Mercury. 
In particular, we present the parameterized post-Newtonian formalism, which has become a useful 
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theoretical framework for analyzing the gravitational experiments within the solar system. This section 
contains a catalog of relativistic gravitational experiments possible with a Mercury Orbiter. We present 
both the quantitative and qualitative analyses of the measurable effects. In section VI, we discuss briefly 
the radio science instrumentation required for the Mercury Or6«ter mission under consideration by ESA. 

2 Mercury Gravity Field 

The determination of Mercury's gravity field could be obtained by the reduction of Doppler radio 
tracking data. This is an old and proven technique used for the moon, Mars and Venus (Konopliv 
and Sjogren, 1994). The field would be modeled with spherical harmonic coefficients does providing a 
smooth global model. Given the power of modern computers, it is feasible to estimate literally thousands 
of harmonic coefficients. The resulting coefficients would be interpreted in terms of a gravity surface 
function, which would be plotted as a surface gravity map. For other planets, in particular Venus, Earth, 
the moon and Mars, the correlation of gravity maps with surface topography and imaging has revealed 
the isostatic state of surface features. For example, similar features at different locations on a planet 
can have significantly different internal density variations or structure. Geophysical interpretation can 
yield information on the distribution of crustal thickness and dynamic hot spots. The degree of isostatic 
compensation may indicate an old crust, nearly completely compensated, or a crust with sufficient 
rigidity to imply a cold interior. 

For Mercury, the second order gravity coefficients C20, ^21, S'21, C22 and S22 would provide the orien- 
tation of the axes of the moments of inertia, which would be compared with the spin axis orientation. A 
more accurate C20 coefficient would provide a better-determined gravitational oblateness, and C22 co- 
efficient would reveal a better-determined gravitational ellipticity of the equator. At a more local level, 
all the other coefficients, to as high a degree and order as possible, would provide a gravity map for 
comparison with surface topography. The most obvious and interesting correlations would be negative 
(i.e., gravity highs in topographic lows or gravity lows in topographic highs). On the moon and Mars 
there are negative correlations for some large basins (but not all). This has placed strong constraints 
on interior thermal models. We anticipate that a gravity map of Mercury's Coloris Basin would place 
significant constraints on its underlying internal structure. 

As a rule of thumb, the spatial resolution of gravity features derived from coherent X-band Doppler 
data would be less than the spacecraft altitude. For a 400 km altitude periherm, the best resolution 
would be approximately 300 km or 7.1° on the surface of Mercury and would correspond to a 25th degree 
and order spherical harmonic model. At a lower 200 km altitude a 50th degree and order model could 
be obtained for a total of 2600 harmonic coefficients. The ideal data acquisition period would be 59 days 
of Doppler data covering at least ±1 hour of periapsis. Data within 5° to 10° of solar conjunction would 
be too noisy for gravity field determination. During gravity passes other non-gravitational forces would 
be minimized. These data would provide complete global coverage. Data over 59 days at a different 
viewing geometry would significantly strengthen the solution. 

In addition to the gravity coefficients, other parameters would be estimated simultaneously from the 
coherent Doppler data. These include the Love number k2, which would provide a measure of whether 
there is a solid or liquid core, translational forces from solar radiation impinging on spacecraft parts. 
Mercury's albedo, corrections to the Mercury ephemeris, and location of the spin pole and rotation rate. 

The final data products that would be produced for geophysical interpretation are: 

1. The spherical harmonic coefficients. 

2. A full covariance matrix on all the coefficients (i.e., 1 — a uncertainties and their 
correlations). 

3. A 1° X 1° digitized vertical surface gravity map. 

4. A 1° X 1° digitized geoid map. 
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5. 1 — (T uncertainty maps on gravity and geoid. 

6. If topography data, 1° x 1° Bouguer gravity maps. 

7. 1° X 1° Isostatic Anomaly map. 

3 Mercury Ionosphere 

Mercury was encountered three times by the spacecraft Mariner 10, but was occulted by the planet's 
disk only at the first flyby on 29 March 1974. The S/X-band dual-frequency radio links from spacecraft 
to Earth were monitored during the flyby (Howard et al, 1974). The attempt to detect the ionosphere 
from the dispersive phase shift between the two downlinks, however, yielded only an upper limit of 
1000 cm~'^ for the electron density in the altitude range from the surface up to 1000 km (Fjeldbo et 
al, 1976). These observations were used to infer an upper limit to the surface density of the Hermean 
atmosphere of 10^ molecules cm~^. Taking a surface temperature of 500 K, this corresponds to a pressure 
of ca. 10~^^ bar (a very decent vacuum!). The mission profile that evolved during the ESA Mercury 
Orbiter assessment study (Hechler, 1994) offers definite prospects for radiometric investigations. The 
geometry has many parallels to that designed for Mars Orbiter, i.e. regularly recurring occupations 
at polar latitudes during certain phases of the mission (Tyler et al, 1992). Whereas the atmosphere 
on Mars is dense enough to be synoptically measured on a global basis. Mercury's high temperature 
and low surface gravity combine to produce only a residual "exosphere" with the above mentioned 
diminutive surface densities. Detecting an ionosphere on Mercury using radio occupation techniques 
is a difficult, but not hopeless, task. Generally speaking, the ionosphere should be easier to detect 
than the neutral atmosphere. One classical example of an ionospheric detection was the Pioneer 10 
measurement during occultation of the Jovian moon lo (Kliore et al., 1974), for which peak densities of 
6 X 10^ electrons cm~^ were derived. The geometry of this particular occultation enabled the detection 
to be made with only an S-band downlink. The ionospheric effect on the S-band downlink during 
the Voyager 2 occultation of Neptune's moon Triton was about five times greater than that of the 
neutral atmosphere (Tyler et al., 1989). There are, of course, counterexamples such as Titan (Lindal 
et al., 1983), which produced dramatic refraction effects from its dense, extended atmosphere, but 
no discernable ionospheric signature. Radio propagation experiments on previous ESA interplanetary 
missions (e.g. Giotto, Ulysses) have utilized the extensive facilities of the NASA Deep Space Network. 
It is expected that similar ground equipment will be available for the Mercury Orbiter. The assessment 
study for the Mercury Orbiter recommended an S /X-Band radio communications system for both up- 
and downlinks (ESA Study Team, 1994). For the dual- frequency (S/X-band) downlink configuration 
on Ulysses, the smallest detectable change in differential phase (i.e. S-band minus 3/11 X-band) was 
estimated to be ca. 5° at S-band, corresponding to an electron content of about 10^^ electrons m~^ 
(Bird et al., 1992). The occultation geometry afforded by an orbiting spacecraft, however, for which 
the lateral motion of the signal ray path through the ionized medium is quite rapid, may still enable 
detection of intervening plasma without the dual-frequency capability. For a geometry of this type, 
the smallest phase change measureable is estimated as ~0.5°. This is also the level of detectability 
attributed to with the uplink radio occultation instrumentation presently being developed for the Pluto 
Express Mission [G.L. Tyler, 1995. private communication]. This tiny change in phase, which provides 
a fair estimate of the present-day detection threshold, corresponds to a neutral column density of 5 x 
10^^ molecules m~^ or an electron columnar content of 10^^ electrons m~^ at X-band. For a spherically 
symmetric Hermean atmosphere (ionosphere), this can be shown to translate to peak number densities 
along the ray path at the planet's limb of 5 x 10^^ molecules cm~^ (100 electrons cm~^). Using the 
upper limits quoted by Fjeldbo et al. (1976), column densities for the ray path tangent to the limb may 
be estimated for Mercury. The neutral column density is found to be 10^^ molecules m~^, which is, 
unfortunately, about six orders of magnitude below the detection threshold. Phase changes induced by 
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the neutral gas are thus probably negligible. The ionized component, however, if we assume the peak 
reaches 10^ electrons cm~^, can produce a content of 10^^ electrons m~^, which would place it above 
the lower limits of detectability. This estimate, it should be emphasized, is only valid for a global-scale 
ionosphere. However, as pointed out during the Mercury Orbiter assessment study (ESA Study Team, 
1994), atomic sodium, which emanates from various hot spots on the planet's surface, may well be the 
most prevalent component of the Hermean exosphere. Mean surface densities of up to 5 x 10*^ cm~^ 
are inferred (Cheng et al, 1987). The Na atom is photoionised on time scales of a few hours, so this is 
presumed to be the most important source of ionospheric plasma. The proximity of the magnetopause 
and the tailward-streaming solar wind, however, may cause many ions to be swept away at rates even 
quicker than those of the sources. If the electron density is concentrated in the atmosphere above a 
local source, densities considerably higher than 1000 cm~^ would be necessary to produce a measureable 
effect. Regularly recurring Earth occultations, during which radio tracking at ground-based stations is 
interrupted, are included in virtually every Mercury Orbiter mission scenario because they are basically 
"unavoidable". Two such occultation intervals, one near inferior conjunction and one near superior 
conjunction, were found to occur for the prototype Mercury Orbiter mission (Hechler, 1994). The total 
number of such occultations, each with ingress and egress opportunities, is well over 100. 

Another feature, as mentioned above, is the rapid swath taken by the radio ray path through the 
ionosphere. This magnifies the plasma Doppler shift on the occulted radio signal, thereby increasing 
the detection probability. Synoptic Doppler observations of ionospheric electron densities on a Mercury 
Orbiter would thus appear to be very difficult, but still worth a good look. Even in the absence of a 
permanent global ionosphere it is still possible, however, that a systematic examination of the radio 
occultation data would result in the serendipitous detection of near-surface plasma clouds. 

4 Mercury Ranging Data 

A summary and discussion of data used in the fundamental ephemcrides DE200/LE200 is available 
(Standish, 1990; Standish et al., 1992). Those data include Mercury radar ranging data spanning the 
years 1966-1974. Analysis of the older data plus additional data spanning the years 1974 to 1990 is 
underway (Anderson et al., 1995). Included are two range fixes from the Mariner 10 Mercury flybys on 
March 29, 1974 and March 16, 1975 (Anderson et al., 1987). The current JPL set of reduced Mercury 
radar ranging data is summarized in Table 1. The Goldstone data are from 34 m and 70 m stations 
located in California's Mojave desert, the Arecibo data are from Puerto Rico, while the Haystack data 
are from Tyngsboro, Massachusetts. 

The data used in current relativity tests arc reduced ranging measurements to Mercury's surface. 
Both radar time delay and Doppler data have been used in the reduction (for an explanation of the 
Doppler-delay technique see e.g. Ingalls and Rainville 1972; Shapiro et al., 1972; Harmon et al, 1986). 
The 1978-1982 reduced data are from archives at the Harvard Smithsonian Center for Astrophysics (J. 
Chandler, private communication). Systematic errors in radar ranging data limit their effectiveness for 
tests of general relativity, most notably the excess precession of Mercury's perihelion, and topography 
models have been introduced for purposes of minimizing the error (Anderson et al., 1991; Pitjeva, 1993; 
Anderson et al, 1995). New data, much less affected by systematic error, will become available with 
Mercury Orbiter. Specifically, transponded ranging, similar to that from the Mariner 9 Mars Orbiter, 
will be available. Eventually, with the introduction of a Mercury Lander, even better surface transponded 
ranging will become available, similar to ranging data from the two Viking Landers. Unfortunately, no 
Venus Orbiter has carried a ranging transponder to date. However, we have removed Venus topography 
from Venus radar ranging data by using the Pioneer 12 radar altimetry measurements, and one range 
fix is available from the Galileo spacecraft flyby in 1990 (Anderson et al., 1991). In addition, the MGN 
X-band Doppler data provide a measure of the line-of-sight velocity to an accuracy of 0.02 mm/s (or 
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~ 10m range accuracy). Of all the inner planets, only Mercury currently requires a parameterized 
topography model for ranging data analysis. The Mercury Orbiter will remove this restriction, and the 
orbital motion of Mercury will be placed on an equal footing with Venus and Mars. 

5 Relativistic Gravitational Tests 

Mercury is the closest to the Sun of all the planets of the terrestrial group and, due to its unique 
location and orbital parameters, it is especially suited to relativistic gravitational experiments. The 
short period of its solar orbit will allow experiments over several orbital revolutions. In this section 
we will define the list of gravitational expermients for the Mercury Orbiter mission. Depending on the 
nature of the measureable effect, there are two main groups of possible experiments, namely: (i) the 
celestial mechanical studies of the orbital motion and rotation of Mercury in the gravitational field of 
the Sun, as well as observations of the motion of the Mercury Orbiter in the solar gravitational field at 
the vicinity of the planet Mercury, and (ii) the studies of the gravitational infiuence on the propagation 
of the radio waves along the Earth-Mercury path. We will discuss separately the possible experiments 
from both these groups and will present the estimates of the magnitude and physical character of the 
measured effects. Analysis, performed in this section, is directed towards the future mission, so we will 
show which relativistic effects can be measured, and how accurately these measurements can be done. 

5.1 Parameterized Post-Newtonian Gravity. 

Metric theories of gravity have a special position among all the other possible theoretical models. The 
reason is that, independently of the many different principles at their foundations, the gravitational field 
in these theories affects the matter directly through the metric tensor of the Riemannian space-time 
9mn, which is determined from the field equations of a particular theory of gravity. As a result, in 
contrast to Newtonian gravity, this tensor contains the properties of a particular gravitational theory as 
well as carries the information about the gravitational field of the bodies. This property of the metric 
tensor enables one to analyze the motion of matter in one or another metric theory of gravity based 
only on the basic principles of modern theoretical physics. 

In order to accumulate the features of many modern metric theories of gravity in one theoretical 
scheme as well as to create the 'versalite instrument' to plan the gravitational experiments and to 
analyze the data obtained, Nordtvedt and Will have proposed the parameterized post-Newtonian (PPN) 
formalism (Nordtvedt 1968a; Will 1971; Will & Nordtvedt 1972). This formalism allows one to describe 
within the common framework the motion of the celestial bodies for a wide class of metric theories of 
gravity. Within the accuracy of modern experimental techniques, the PPN formalism provides a useful 
foundation for testing the predictions of different metric theories of gravity within the weak gravitational 
field and slow motion approximation appropriate for the solar system. 

The metric tensor of the general Riemannian space-time in the PPN formalism is generated by some 
given distribution of matter in the form of an ideal fiuid. It is represented by a sum of gravitational 
potentials with arbitrary coefficients, the PPN parameters. The most general form of this metric tensor 
in four dimensions may be written as followffl: 

500 = 1 - 2C/ + 2{/3 + t)[/2 - (27 + 2 + as + Ci)"^! + Ci^ + 2(C + r)$^- 



(37 + 1 - 2/3 + r + C2)^2 + (1 + C3)^3 + 3(7 + C4)^4 + i^X, 



00 



^The geometrical units Ti = c = G = 1 are used throughout as is the metric convention (+ 
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-(ai - 02 - a3)wi^w^'U + a2w''w''U^^ + (ai - 2a3)w''Vf, + ©(c"*^), (la) 
5oa = ^(47 + 3 + 2z/ + 2r + ai - a2 + Ci)K.+ 
+ ^(1 - 2i/ - 2r + a2 - Ci)W^a - ^(ai - 2a2)w^aC/ + a2wt'Ua„ + 0(c"^), (16) 

9aP = 7a/3(l + 2(7 + r)U) + 2TUaf} + 0{c~^). (Ic) 
with the generalized gravitational potentials given as follows: 

J \Z Z- \ J \z z \ 



U^\z^) = I poiz'n ^'^ "'^^ d'z"^, 



J \^ ^ \ J \z z \ 

J \z Z \ J \z z \ 

^U^n />0(-'^)P0(."-)^^f^[||-^ - ^^^]d3.-d3.-. (2) 

Besides the PPN parameters (7, /3, C, 01—03; Ci— C4)) the expression (1) contains two other parameters 
v and r (Denisov and Turyshev, 1990; Turyshev, 1990). The parameter v reflects the specific choice 
of the gauge conditions and for the standard PPN gauge it is given as u = 0, but for harmonic gauge 
conditions one should choose v = \- The parameter r describes the possible anisotropy of the space-time 
and it corresponds to a different spatial coordinates, which may be chosen for modelling the experimental 
situation. For example, the case r = corresponds to harmonic coordinates, while r = — 1 corresponds 
to the standard (Schwarzschild) coordinates. Then, a particular metric theory of gravity in the PPN 
formalism with a specific coordinate gauge (z^, r) might be fully characterized by the means of the ten 
PPN parameters. This formalism uniquely prescribes the values of these parameters for each particular 
theory under study. In the standard PPN gauge (i.e. in the case when 1/ = r = 0) these parameters 
have clear physical meaning. The parameter 7 represents the measure of the curvature of the space- 
time created by the unit rest mass; the parameter /? is the measure of the non-linearity of the law of 
superposition of the gravitational fields in a theory of gravity (or the measure of the metricity); the 
parameter C, quantifies the violation of the local position invariancy; the group of parameters 01,02,0(3 
specify the violation of Lorentz invariancy (or the presence of the privileged reference frame), and, finally, 
the parameters Ci) C2) Cs; C4 reflect the violation of the law of total momentum conservation for a closed 
gravitating system. Note that general relativity, when analyzed in standard PPN gauge, gives: 7 = /3 = 1 
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and all the other eight parameters vanish. Whereas for the Brans-Dicke theory, (3 = 1,7 = where 
ui is an unknown dimensionless parameter of this theory. 

The main properties of the solution (l)-(2) are well established and widely in use in modern astro- 
nomical practice (Moyer, 1971; Moyer, 1981; Brimibcrg, 1991; Standish et al, 1992; Wih, 1993). For 
practical purposes one chooses the inertial reference frame located in the center of mass of an isolated 
distribution of matter, then, by performing a power expansion of the potentials in terms of spherical 
harmonics, one obtains the post-Newtonian definitions of the mass of the body nia, it's center of mass 
Z^, momentum and angular momentum S*"^. Thus the definitions for the mass nia and coordinates 
of the center of mass of the body Z" in any inertial reference frame are given by the formulae (for more 
detailed analysis see (Damour, 1983) and (Will, 1993) and references therein): 

m„ = ^ dh': p*o (l + n + if/ - ^vxv^ + 0(c-^)) , (3a) 

Ula J a ^ Z Z ' 

where the conserved mass density is defined by /9*g = paQ\/—gu^- 

It has been shown (Fock, 1959; Ni & Zimmerman, 1978) that for an isolated distribution of matter in 
this approximation there exist a set of inertial reference frames and ten integrals of motion corresponding 
to ten conservation laws. Moreover, Ni and Zimmerman have shown (1978) that in order for a metric 
theory to have the total energy- momentum conservation laws (so-called fully conservative theories), 
only three PPN parameters may have non-zero means, namely 7,/3 and ^. For the future analysis of 
gravitational experiments on Mercury Orbiter, we will limit ourselves to this class of metric theories of 
gravity. 

In the post-Newtonian approximation the mass rua (3a) is conserved and the centre of mass Z^ moves 
in space with a constant velocity along a straight line. Thus, Z^{t) = ■ t + /C", where the constants 

= dZ^/dt and /C° are the body's momentum and center of inertia, respectively. One can choose 
from the set of inertial reference frame the barycentric inertial one. In this frame the functions Z"^ must 
equal zero for any moment of time. This condition can be satisfied by applying to the metric eqs.(l) the 
post-Galilean transformations (Chandrasekhar &: Contopulos, 1967) where the constant velocity and 
displacement of origin 6° should be selected in a such a way that "P" and /C" equal zero (for details see 
Kopejkin, 1988; Will, 1993). The resulting barycentric inertial reference frame has been adopted for the 
fundamental planetary and lunar ephemeridcs (Ncwhall et al, 1983; Standish et al., 1992). Moreover, 
the coordinate time of the solar barycentric (harmonic) reference frame is the TDB time scale, which 
has been adopted in modern astronomical practice (Fukushima, 1995). 

In order to analyze the motion of bodies in the solar system barycentric reference frame one may 
obtain the Lagrangian function describing the motion of N extended sclf-graviting bodied (Tury- 
shev, 1990; Will, 1993). Within the accuracy necessary for our future discussion of the gravitational 
experiments on the Mercury Orbiter mission, this function can be presented in the following form: 

^ m 1 

a 

N N /, 

-EE ^ + (3 + 7 - 4/3 + 3C)J)„ - CnMl,.^'a'- 

a b^a "^"-^ 

1 3 

-{l + r+ -)Va,,Va^ + (7 + T + -)Va^Vb^'- 
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Tab 



1 1 , /V\ 

(7 + - (7 + ^Hi, + nabxriab^,^ + 

^ """ab J 



^ n 'ab' 



N N N 



bi^a 



N 



^ab "ab^ac-^ „ 



(4) 



o- b^a c^a,b 

where ma is the isolated rest mass of a body a, the vector r" is the barycentric radius-vector of this 
body, the vector r"^ = — is the vector directed from body a to body b, and the vector n^*^ = r^ii/vab 
is the usual notation for the unit vector along this direction. It should be noted that the expression (4) 
does not depend on the parameter v, which confirms that this parameter is the gauge parameter only. 
The tensor /^"^ is the spatial trace-free (STF) (Thorne, 1980) tensor of the quadrupole moment of body 
a defined by 

ir = 2^ / d'^'Ploi^iK'^' - r'^^Wa^) = - r''Ja2, (5) 

where the quantity Ja2 is the quadrupole coefficient. The tensor Sf^" is the body's intrinsic STF spin 
moment which is given as: 

= ^ / d'^aPloMK^a - (6) 
'''■a Ja 

where -u^ is the rotational velocity of body a. Finally, the tensor Jl^'^ is the body's gravitational binding 
energy tensor: 



a/3 



2ma 



I -x/fc" 1 3 



(7) 



Furthermore, with the help of the Lagrangian (4) , one can obtain the corresponding equations of motion 
of an arbitrary extended body a as follows: 



+ E 



nib 



67^0 a& 67^0 ab 



'Kb + — + ^- 

"""ab 



n 



ab 

rab 



{2/3 + 27 + 2r + l)ma + (2/3 + 27 + 2T)mb) 
+ mbm,V2b, + 0{c-% 



+ 



(8) 



where, in order to account for the influence of the gravitational binding energy ^'^^ , we have introduced 
in the equation (8) the tensor of passive gravitational rest mass Mf as follows 



M 



a/3 



rub 



j«/3 (1 + (3 + ^ _ + ^Q^^ _ ^(^nabxnabt^^l^) + K^f + 0{c-^) 



(9) 
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and the unit vector has also been corrected by the gravitational binding energy and the tensor of 
the quadrupole moment I^^ of the body a under question: 



(1 
1 + (3 + 7 - 4^ + SQ^a - Knabxriab^^a'' + 5nabxnab^f2 



ab ^ 



ab 



The term yl"^ in the expression (8) is the orbital term which is given as follows: 

= v>abx{v^ - (27 + 2r + 1)^;^,) + 



+Kb{vaXVa - (7 + 1 - T')VabxVab + ^^inabx^ab^f " ^^riabxVbf), 

The spin-orbital term has the form: 



B2b = il + 27)vabx{S^^ + S^^) + \vax{S^' - S^^)+ 



(10) 



(11) 



+-(1 + 2j)nabxvabf, [niiS^^ + S^^) - n„;."(5f + ) 



+ 



3 

+ 2^abX 



KbiVa^S^^ - VbpS^^) + nabpvJS\ 
The term is caused by the oblateness of the bodies in the system: 



•aX 



Cab = '^riabfjlb'^ + ^Kb^abX^abult^ ■ 



(12) 



(13) 



And, finally, the contribution D'^^^^ to the equations of motion (8) of body a, caused by the interaction 
of the other planets {h ^ a,c ^ a, h) with each other is: 



abc 



ab 
''ab 



(1-2/3 + 20— -2(/? + 7-C) — 

^6c "^ac 



-(C + r)^inbcx + ncax) + (C - r)^^+ 



' ab 



+ l(l + 2C-2r)^ 



^ + 2(1+7)^ 



(14) 

'be ' ''^bc'^a.b 

where the tensors A(^^ = jf^" + n'^n'^ and 11^^ = + 3n^^n^(, are the projecting and polarizing 
operators respectively. 

The metric tensor (1), the Lagrangian function (4) and the equations of motion (8) define the behavior 
of the celestial bodies in the post-Newtonian approximation. Hence, they may be used for producing the 
numerical codes in relativistic orbit determination formalisms for planets and satellites (Moyer, 1981; 
Ries et al, 1991; Standish et al, 1992) as well as for analysing the gravitational experiments in the 
solar system (Will, 1993; Pitjeva, 1993; Anderson et al, 1996). 
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5.2 Metric Tests of Parametrized Gravitation. 

By means of a topographic Legendre expansion complete through the second degree and order, the 
systematic error in Mercury radar ranging has been reduced significantly (Anderson et al, 1995). How- 
ever, a Mercury Orbiter is required before significant improvements in relativity tests become possible. 
Currently, the precession rate of Mercury's perihelion, in excess of the 530 arcsec per century {"/cy) 
from planetary perturbations, is 43.13 " /cy with a realistic standard error of 0.14 " /cy (Anderson et 
al., 1991). After taking into account a small excess precession from the solar oblateness, Anderson et 
al. find that this result is consistent with general relativity. Pitjeva (1993) has obtained a similar result 
but with a smaller estimated error of 0.052 "/cy. Similarly, attempts to detect a time variation in the 
gravitational constant G using Mercury's orbital motion have been unsuccessful, again consistent with 
general relativity. The current result (Pitjeva, 1993) is G/G = (4.7 ±4.7) x 10"^^ yr"^ 

Metric tests utilizing a Mercury Orbiter have been studied both at JPL and at the Joint Institute 
for Laboratory Astrophysics (JILA) and the University of Colorado. The JPL studies, conducted in the 
1970's, assumed that orbiter tracking could provide daily measurements (normal points) between the 
Earth and Mercury centers of mass with a 10 m standard error. A covariance analysis was performed 
utilizing a 16-parameter model consisting of six orbital elements for Mercury and Earth respectively, 
the metric relativity parameters /? and 7, the solar quadrupole moment J20, and the conversion factor 
A between unit distance (Astronomical Unit) and the distance in meters between Earth and Mercury. 
It was assumed that no other systematic effects were present, and that the normal-point residuals after 
removal of the 16-parameter model would be white and Gaussian. The total data interval, assumed 
equal to two years, corresponded to 730 measurements. Under the assumed random distribution of data, 
the error on the mean Earth-Mercury distance was 10/\/730 = 37 cm. The JPL studies showed, based 
on a covariance analysis, that the primary metric relativity result from a Mercury Orbiter mission would 
be the determination of the parameter 7, the parameter that measures the amount of spatial curvature 
caused by solar gravitation. The standard error was 0.0006, about a factor of two improvement over the 
Viking Lander determination. This accuracy reflected the effect of spatial curvature on the propagation 
of the ranging signal and also its effect on Mercury's orbit, in particular the precession of the perihelion. 
The error in the metric parameter /3 and the error in the solar J20 were competitive with current results, 
but not significantly better. 

Within the last five years, a more detailed covariance analysis by the JILA group (Ashby et al., 
1995) assumed 6 cm ranging accuracy over a data interval of two years, but with only 40 independent 
measurements of range. Unmodeled systematic errors were accounted for with a modified worst-case 
error analysis. Even so, the JILA group concluded that a two-order of magnitude improvement was 
possible in the perihelion advance, the relativistic time delay, and a possible time variation in the 
gravitational constant G as measured in atomic units. However, the particular orbit proposed by ESA 
for its 2000 Plus mission was not analyzed. It is almost certain that the potential of the ESA mission 
lies somewhere between the rather pessimistic JPL error analysis and the JILA analysis of an orbiter 
mission more nearly optimized for relativity testing. 

5.2.1 Mercury's Perihelion Advance. 

The secular trend of Mercury's perihelion in the four-dimensional (4D) PPN formalism depends on the 
linear combination of the PPN parameters 7 and /3 and the solar quadrupole coefficient J20 (Nobili and 
Will, 1986; Heimberger et al, 1990; Will, 1993): 

7r4D = (2 + 27-^)^ 2-T + 7( In cos im - 1), cy (15a) 

' c^aMil- eir) 4\aM^ (1 - el.)^ 
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where um, a.nd bm are the mean distance, mean motion, inchnation and eccentricity of Mercury's 

orbit. The parameters Mq and Rq are the solar mass and radius respectievely. For the Mercury's orbital 
parameters one obtains: 

TTiD = 42'(98 \\{2 + 2j- P) + 0.296 • J20 x 10^1 , "/cy (156) 

Thus, the accuracy of the relativity tests on the Mercury Orbiter mission will depend on our knowl- 
edge of the solar gravity field. The major source of uncertainty in these measurements is the solar 
quadrupole moment J2q - As evidenced by the oblatcncss of the photosphere (Brown et al, 1989) and 
perturbations in frequencies of solar oscillations, the internal structure of the Sun is slightly asphcrical. 
The amount of this asphericity is uncertain. It has been suggested that it could be significantly larger 
than calculated for a simply rotating star, and that the internal rotation rate varies with the solar cycle 
(Goode and Dziembowski 1991). Solar oscillation data suggest that most of the Sun rotates slightly 
slower than the surface except possibly for a more rapidly rotating core (Duvall and Harvey, 1984). 
An independent measurement of J20 performed with the Mercury Orbiter would provide a valuable 
direct confirmation of the indirect helioseismology value Furthermore, there are suggestions of a rapidly 
rotating core, but helioseismology determinations have been limited by uncertainties for depths below 
0.4 solar radii (Libbrecht and Woodard, 1991). 

The Mercury Orbiter will help us understand this asphericity and independently will enable us to 
gain some important data on the properties of the solar interior and the features of it's rotational 
motion. Preliminary analysis of a Mercury Orbiter mission suggests that J20 can be determined to at 
best ~ 10~^ (Ashby, Bender and Wahr, 1995). Even a determination ten times worse, on the order of 
10% of the total effect, would be comparable in accuracy to the solar oscillation determination (Brown 
et al, 1989) 

In addition to this, the studies of Mercury's perihelion advance might provide us with an opportunity 
to test the hypothesis of the multi-dimensionality (D > 4) of our physical world (Gross and Perry, 1983; 
Lim and Wesson, 1992; Kalligas et al., 1995). As it was shown by Lim and Wesson, multi-dimensional 
extensions of general relativity may produce observable contributions to the four-dimensional physical 
picture. It has been conjectured that even a small extra (curved) dimension may affect the dominant 
terms in the four-dimensional metric, thereby altering the theoretical value of the perihelion shift. 
For example, in the five-dimentional Kaluza-Klein theory, the predicted value of Mercury's perihelion 
advance is given as 

TTSD = 7r4D(l + k^), (15c) 

where k is the small parameter representing the motion of the body in the extra fifth dimension. 

It should be noted that the Mercury Orbiter itself, being placed on orbit around Mercury, will 
experience the phenomenon of periapse advance as well. However, we expect that uncertainties in 
Mercury's gravity field will mask the relativistic precession, at least at a level of interest for ruling out 
alternative gravitational theories. 

5.2.2 The Precession Phenomena. 

In addition to the perihelion advance metric theories predict several precession phenomena associated 
with the angular momentum of the bodies. Thus, according to PPN formalism, the spin of a gyroscope 
So attached to a test body orbiting around the Sun processes with respect to a distant standard of rest 
- quasars or distant galaxies. The theory of Fermi- Walker transport describes the motion of the spin 
vector of a gyroscope by the relation 

^ = [nxs-h]. (16) 
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where the angular velocity Ci has the following form (Will, 1993): 

n = -l[vxa] + ij + ^)[v X VUq]- 

-(1 +7)0(5; -3(S;n)n). (17) 

Here v is the velocity of the body, a = dv/dr — VUq is the non-gravitational acceleration of the body 
under consideration (for example from the solar radiation pressure) , r the body's heliocentric radius 
vector, and Uq is the gravitational potential of the Sun at the body's location. 

The first term in (17) is the Thomas precession. It is a special relativistic effect due to the non- 
comutativity of nonaligned Lorentz transformations (Thomas, 1927; Bini et al, 1994). It also can be 
interpreted as a coupling between the bodies velocity and the non-gravitational forces acting on it. The 
studies of this precession may provide an additional knowledge of the non-gravitational environment of 
Mercury and its interaction with the solar wind. As for the quasi-Lorcntz transformations to the proper 
coordinate reference frame, this term produces a negligible small contribution to the total precession 
rate of the satellite's orbit. 

The second term is known as a geodetic precession (De-Sitter, 1916). This term arises in any non- 
homogeneous gravitational field because of the parallel transport of a direction defined by sq. It can be 
viewed as spin precession caused by a coupling between the particle velocity v and the static part of the 
space-time geometry (1). For Mercury orbiting the Sun this precession has the form: 

^G = {l+l)^{fxv), (18) 

where /Iq is the solar gravitational constant, r is the distance from Mercury to the Sun and v is its 
the orbital velocity. This effect should be studied for the Mercury Orbiter, which, being placed in orbit 

around Mercury is in effect a gyroscope orbiting the Sun. Thus, if we introduce the angular momentum 
per unit mass, L = f x v, oi Mercury in solar orbit, the equation (18) shows that da is directed along 
the pole of ecliptic, in the direction of L. The vector has constant part 

Jlo = J(l + 27)^ = i±^- 0.205 7yr, (19a) 

with the significant correction due to the eccentricity cm of the Mercury's orbit, 

Hi cosa;M^ = ^(1 + 27) ^^^^ Cm cos UMt = ^ " 0.126 cos cjm^ "/y^, (19^) 
2 qm 3 

where ujm is Mercury's siderial frequency and t is reckoned from a perihelion passage; om is the semi- 
major axis of Mercury's orbit. 

This effect has been studied for the motion of lunar perigee and the existence of the geodetic pre- 
cession was confirmed with an accuracy of 10% (Bertotti, Ciufolini and Bender, 1987). Two other 
groups has analyzed the lunar laser-ranging data to estimate the deviation of the lunar orbit from the 
predictions of general relativity (Shapiro et al, 1988; Dickey et al., 1989). These predictions have been 
confirmed within the standard deviation of 2%. 

Certainly this experiment, being performed in the vicinity of the Mercury, requires complete analysis, 
taking into account the real orbital parameters of the spacecraft. Our preliminary analysis shows that 
the precession of the satellites's orbital plane should include a contribution of order 0.205 "/yr from the 
geodetic precession. Therefore special studies of this precession should be included in future studies of 
the Mercury Orbiter mission. 
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The third term in the expression (17) is known as Lense-Thirring precession J^lt- This term gives 
the relativistic precession of the gyroscope's spin sq caused by the intrinsic angular momentum S of the 
central body. This effect responsible for a small perturbation of orbits of artificial sateUites around the 
Earth (Tapley et ai, 1972; Ries et ai, 1991). Unfortunately, preliminary studies have shown that this 
effect is also negligibly small for the satellite's orbit around Mercury and, moreover, it is masked by the 
uncertainties of the orbit's inclination. 



5.2.3 The Test of the Time Dependence of the Gravitational Constant. 

As pointed out by Dirac, the age of the universe, Hq^, in terms of the atomic unit of time, e^/mc^, is 
of order 10^^, as is the ratio, / (Gnipm) , of the electrical force between the electron and proton to the 
gravitational force between the same particles. This, according to Dirac, suggested that both ratios are 
functions of the age of the universe, and that the gravitational constant, G, might vary with time. In 
order to account for this possibility, and assuming that G = dG/dt is constant in first approximation, 
we define a parameter \q as follows: G/G = —XgHq. Then, the effect of \q on the motion of bodies in 
the solar system can be introduced as an additional perturbation in the barycentric equations of motion 
(8) in the following form: 

5xa^= XgHo ^fiab + 0{t'). (20) 

The major contribution to this effect for the motion of celestial bodies in the solar system is the Sun. 
The importance of the Mercury Orbiter mission will be in providing the tracking data necessary for 
establishing a more accurate determination of the orbital elements of this planet. Then, combined 
together with other solar system data (Chandler et al, 1994; Anderson et ai, 1996), this information 
will enable one to perform a more accurate test of the hypothesis (20) . 



5.2.4 The Planetary Test of the Equivalence Principle. 

The development of the parameterized post-Newtonian (PPN) formalism, has provided a useful frame- 
work for testing the violation of the Strong Equivalence Principle (SEP) for extended bodies (Anderson 
et al, 1996). In that formalism, the ratio of passive gravitational mass mb(g) to inertial mass mb(j) is 
given by 

^ = 1 + V^, (21) 

where mi, is the rest mass of a body h and Q,}, = ^^u^'^^ is its gravitational binding energy (2.7). 
Numerical evaluation of the integral of expression (7) for the standard solar model (Ulrich, 1982) gives 

^ ^ 3.52 • 10~^ (22) 



The SEP violation is quantified in eq.(21) by the parameter 77. In fully-conservative, Lorentz-invariant 
theories of gravity the SEP parameter is related to the PPN parameters by 

r; = 4/3 - 7 - 3 (23o) 

and is more generally related to the complete set of PPN parameters through the relation (Will, 1993) 

r/ = 4;3-7-3-ye-ai + ^ (2^2 - 2Ci - C2). (236) 

A difference between gravitational and inertial masses produces observable perturbations in the 
motion of celestial bodies in the solar system. By analyzing the effect of a non-zero 77 on the dynamics 
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of the Earth-Moon system moving in the gravitational field of the Sun, Nordtvedt (1968b) found a 
polarization of the Moon's orbit in the direction of the Sun with amplitude 6r ~ rjCo, where Cq is a 
constant of order 13m (Nordtvedt effect). The most accurate test of this effect is presently provided by 
Lunar Laser Ranging (LLR) (Shapiro et al, 1976; Dickey et al, 1989), and in the most recent results 
(Dickey et al, 1994; Williams et al. 1995) the parameter rj was determined to be 

ri = -0.0005 ± 0.0011 . (24) 

Also results are available from numerical experiments with combined processing of LLR, spacecraft 
tracking, planetary radar and Very Long Baseline Interferometer (VLBI) data (Chandler et al. 1994). 
Recently, the analysis of the accuracy for planetary SEP violation in Sun-Jupiter-Mars-Earth system 
has been done for future Mars missions (Anderson et al., 1996). We would like to emphasize that 
a measurement of the sun's gravitational to inertial mass ratio can be performed using the Mercury 
ranging experiment. Indeed, let us consider, for example, the dynamics of the four-body Sun-Mercury- 
Earth-Jupiter system in the solar system barycentric inertial reference frame. The quasi-Newtonian 
acceleration of Mercury (M) with respect to the Sun (5), fsM, is straightforwardly calculated from the 
equations (8) to be: 



rsM = rM-rs = -rngM • -3— + b~ " ^ ^"^V:;^ )s^^^' 



rjs rjM 



JS ' JE 



^ \ rjs 



'JS 



' SM 

where m*g}^^ = ms + tum + f] msi^:^^^ + m^^;^)^ • The first and second terms on the right side of 
the equation (25) are the Newtonian and the tidal acceleration terms respectively. Wc will denote the 
last term in this equation as A^. This is the SEP acceleration term which is of order c^^ and it is treated 
as a perturbation on the restricted three-body problem. While it is not the only term of that order, 
the other post-Newtonian terms (suppressed in eq.(25)) do not affect the determination of 77 until 
the second post-Newtonian order (~ c~^). The corresponding SEP effect is evaluated as an alteration 
of the planetary Keplerian orbit. The subscripts {E) and (J) indicate Earth and Jupiter, respectively. 
To good approximation the SEP acceleration A^^ has constant magnitude and points in the direction 
from Jupiter to the Sun, and because it depends only on the mass distribution in the Sun, the Earth 
and Mercury experience the same perturbing acceleration. The responses of the trajectories of each of 
these planets due to the term A^^ determines the perturbation in the Earth-Mercury range and allows a 
detection of the SEP parameter r/ through a ranging experiment. 

The presence of the acceleration term A^j in the equations of motion (25) results in a polarization 
of the orbits of Earth and Mercury, exemplifying the planetary SEP effect. Wc have examined this 
effect on the orbit of Mercury by carrying out first-order perturbation theory about the zeroth order 
of Mercury's orbit, taken to be circular. Then the SEP perturbation produces the polarization of the 
Mercury's orbit of order: Sr = tjCm, where Cm is a constant of order 81 m. However, it turns out that 
the eccentricity correction played a significant role in the similar problem of Mars motion (Anderson 
et al., 1996). One reason for this is that these corrections include "secular" matrix elements which arc 
proportional to the time t. Such elements dominate at large times, and the eccentricity corrections 
thereby qualitatively change the nature of the solution in the linear approximation. Thus, for a given 
value of SEP parameter 77, the polarization effects on the Mercury orbit are more than two orders of 
magnitude larger than on the lunar orbit. The additional result of these studies is that the mass of 
Jupiter, mj, can be determined more accurately from a few years of Earth- Mercury ranging than from 
Pioneer 10,11 and Voyager 1,2 combined. This analysis shows a rich opportunity for obtaining new 
scientific results from the program of ranging measurements to Mercury during the Mercury Orbiter 
mission. Certainly, for the future analysis of the general planetary SEP violation problem one should 
modify the theoretical model to include effects due to Venus, Mars and Saturn and perform numerical 
experiments with combined data collected from the planetary missions, LLR and VLBL 
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Moreover, the analysis of Mercury Orbiter ranging data might provide the opportunity for another 
fundamental test, namely a solar system search for dark matter. The data obtained during this mission 
will allow one to study the conditions leading to the violation of the strong equivalence principle in the 
motion of extended bodies in the solar system resulting from possible composition, shape and rotation 
dependent coupling of dark matter to the matter of different astrophysical objects. 

5.2.5 The Redshift Experiment. 

Another important experiment that could be performed on a Mercury Orbiter mission is a test of the solar 
gravitational redshift, provided a stable frequency standard is flown on the spacecraft. The experiment 
would provide a fundamental test of the theory of general relativity and the Equivalence Principle upon 
which it and other metric theories of gravity are based (Will, 1993). Because in general relativity the 
gravitational redshift of an oscillator or clock depends upon its proximity to a massive body (or more 
precisely the size of the Newtonian potential at its location), a frequency standard at the location of 
Mercury would experience a large, measurable redshift due to the Sun. Moreover, the eccentricity of 
Mercury's orbit would be highly effective in varying the solar potential at the clock, thereby producing 
a distinguishing signature in the redshift. The anticipated frequency variation between perihelion and 
aphelion is to first-order in the eccentricity: 



This contribution is quite considerable and is calculated to be {Sf / fo)eM — 1-1 ^ 10~^. It's absolute 
magnitude, for instance, at the wave length Aq = 3 cm (/o = 10 GHz) is {Sf)ej^ = 110 Hz. We would 
also benefit from the short orbital period of Mercury, which would permit the redshift signature of the 
Sun to be measured several times over the duration of the mission. Depending upon the stability of the 
frequency standard, a Mercury Orbiter redshift experiment could provide a substantial improvement 
over previous tests of the solar redshift. 

At this time, there exist two comparable tests of the solar redshift. The first test was performed 
by using the well-known technique of measuring shifts in the positions of spectral lines of atoms in the 
solar atmosphere. Recent measurements of an oxygen triplet have verified the solar redshift predicted 
by general relativity to an accuracy of 1% (LoPresto et al, 1991). The first test of the solar redshift 
with an interplanetary probe was performed recently with the Galileo spacecraft during its mission to 
Jupiter (Krisher, 1993). In order to obtain gravity-assists from the Earth and Venus, the spacecraft 
traveled on a trajectory that took it in and out of the solar potential. By measuring the frequency of 
an ultra-stable crystal oscillator (USO) flown on the Galileo spacecraft during this phase of the mission, 
the redshift prediction was again verified to an accuracy of 1%. The accuracy obtained was limited 
by the long-term stability of the USO. After calibration of systematic errors (e.g., a linear drift due to 
aging), the fractional frequency stability was roughly 10~^^. 

During a Mercury Orbiter mission, the variation in the solar potential between perihelion and aphe- 
lion would produce a redshift of roughly 1 part in 10^ according to general relativity (for a detailed 
analysis of the anticipated redshift, see Krisher, 1993). With a fractional frequency stability of 10^^^, 
we could test the redshift with an accuracy of 1 part in 10^ (or 0.01%). This level of accuracy is com- 
parable to that obtained by the 1976 Vessot/NASA Gravity-Probe A (GP-A) experiment performed 
in the gravitational field of the Earth with a hydrogen maser on a Scout rocket (Vessot et al., 1980). 
Greater sensitivity could be obtained if an atomic frequency standard were flown on Mercury Orbiter. 
An atomic standard having a typical stability of 1 part in 10^^ would permit the redshift to be tested 
to 1 part in 10^. 




(26) 
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5.2.6 The Shapiro Time Delay. 



According to general relativity, an electromagnetic signal propagating near the Sun will be delayed by 
solar gravity (Shapiro, 1964). In PPN formalism this additional coordinate time delay, Ar, is given by 



where ve and Vp are the respective distances of the Earth and the planet from the Sun, and VEp = Tp — rE 
is the earth-planet distance. At superior conjunction this formula can be approximated by 



where d is the closest distance from the center of the Sun to the beam trajectory. 

By ranging to Mariner, Viking, and Voyager spacecraft during solar conjunction, it has been possible 
to test this effect to a few percent or better (Krisher et al., 1991 and references therein). From maximum 
elongation to superior conjunction, Ar increases from 15 to 240 microseconds for Mercury. However, all 
previous spacecraft tests have depended on an S-band uplink. The Mercury Orbiter may provide the first 
opportunity for an X-band test. The higher frequency, by the factor 11/3, will reduce the systematic 
error from propagation of the radio beam through the solar corona. It will also be possible to perform a 
Doppler relativity test as suggested by Bertotti and Giampieri (1992). The multi conjunctions may be 
useful in constraining certain alternative theories of gravitation involving anisotropic metric potentials. 
The Doppler shift caused by solar gravity could depend on the geometry of the conjunction. 

6 Radio Science Instrumentation 

The Mercury Orbiter radio science objectives require on board instrumentation with the following ca- 
pabilities: 

1. Up-link at X-band (~ 3.5 cm wavelength) from ground stations to the spacecraft. 

2. X-band down-link coherently referenced to the up-link signal. 

3. Crystal or atomic ultra-stable oscillator (USO). 

4. Simultaneous down-link at S (~ 13 cm) and X-bands referenced to the USO. 

Any S-band up-link capability, perhaps provided for commanding or emergency communications, 
would not be used during radio science observations. The S-band up-link provides far superior per- 
formance. Up-link and down-link at higher frequencies (i.g. KA-band) should be considered for the 
Mercury Orbiter mission for purposes of enhancing the gravity and relativity measurements. A multi- 
link Doppler systems should be considered as well. Oven controlled USO's has been included in several 
radio science systems in the past Including Voyager, Galileo and Cassini. The Mercury Orbiter mission 
is far enough in the future that the crystal oscillator will be replaced by an atomic standard. 

Starting with an X-band system, the cost of adding a dual-frequency downlink capability to enhance 
radio propagation measurements would require about 15-20 watts of additional power and 4-5 kg of 
additional mass for an S-band transmitter. Observations could then be performed at both occultation 
ingress and egress by transmitting both frequencies coherently from a common on-board oscillator. It 
would not be necessary to add a special ultra-stable oscillator (USO) for this purpose as long as the 
downlinks are phase coherent. 



(Will, 1993; Lebach et al, 1995) 





(27) 
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Appendix A: The astrophysical parameters used in the paper. 

Solar radius : Rq = 695, 980 km, 

Solar gravitational constant : /x© = ^ ® = 1.4766 km, 

Solar quadrupole coefficient (Brown et ai, 1989) : J20 = (1.7 ± 0.17) x 10"''', 
Solar rotation period : r© = 25.36 days. 
Mercury's mean distance : qm = 0.3870984 AU = 57.91 x 10^ km. 
Mercury's radius : Rm = 2, 439 km, 

Mercury's gravitational constant : /xm = ^^J^ = 1.695 x 1O~^/X0, 
Mercury's sidereal period : Tm = 0.241 yr = 87.96 days. 
Mercury's rotational period : tm = 59.7 days , 
Eccentricity of Mercury's orbit : cm = 0.20561421, 
Jupiter's gravitational constant : fij = 9.547 x lO~^/i0, 
Jupiter's sidereal period : Tj = 11.865 yr, 
Astronomical Unit : AU = 1.49597892(1) x 10^^ cm. 
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Mercury Radar-Ranging 


; Measurements 


Timespan 


Antenna 


Number 






of Observations 


1967-1971 


Arecibo 


85 


1966-1971 


Haystack 


217 


1971-1974 


Goldstone 


38 


1974-1975 


Mariner 10 


2 


1978-1982 


Arecibo 


157 


1986-1990 


Goldstone 


132 
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